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Abstract. A recently published study on long term evolution of the frequencies of the kilohertz quasi- 
periodic oscillations (QPOs) in the atoll source 4U 1636-53 concluded that there is no preferred frequency 
ratio in a distribution of twin QPOs that was inferred from the distribution of a single frequency alone. How- 
ever, we find that the distribution of the ratio of actually observed pairs of kHz QPO frequencies is peaked 
close to the 3/2 value, and possibly also close to the 5/4 ratio. To resolve the apparent contradiction between 
the two studies, we examine in detail the frequency distributions of the lower kHz QPO and the upper kHz 
QPO detected in our data set. We demonstrate that for each of the two kHz QPOs (the lower or the upper), 
the frequency distribution in all detections of a QPO differs from the distribution of frequency of the same 
QPO in the subset of observations where both the kHz QPOs are detected. We conclude that detections of 
individual QPOs alone should not be used for calculation of the distribution of the frequency ratios. 
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1 Introduction 

Kluzniak & Abramowicz (2000) suggested that the kHz twin peak QPOs, observed in the Fourier power 
spectra (PDS) from accreting neutron stars, originate in a non-linear resonance that is possible only in strong 
gravity. It was reported later (Abramowicz et al. 2003) that the ratio v L /v v of the upper and lower QPO 
frequency in neutron stars usually clusters close to the rational ratio 2/3, with some frequency pairs possibly 
clustering close to other ratios, such as 0.78. 

Belloni et al. (2005) re-examined the study of Abramowicz et al. (2003) for a larger set of detections of 
a single kHz QPO and, on the assumption of a correlation between the observed QPO and the unobserved 
second QPO, confirmed that their (inverse) frequency ratio Vu/v L would cluster most often close to the 3/2 
value and less often close to other rational numbers (e.g., 5/4 and 4/3). Because a distribution of the ratios of 
two correlated quantities is largely determined by the distribution of either one of them, Belloni et al. (2005) 
argued that the peaks in the distribution of kHz frequency ratios reported by Abramowicz et al. (2003) reflect 
peaks of unknown origin in the distribution of a single (upper or lower) kHz QPO. Further, they argued that 
such clustering does not provide any useful information about a possible underlying physical mechanism. 
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Fi gure 1 : Lett: Frequency histogram for the lower QPO. The subset of lower QPOs detected simultane- 
ously with the upper QPO is denoted by darker bars, which are labeled "lower in twin". Right Analogous 
histograms for the upper QPO. 



A more recent study of Belloni et al. (2007) is based on a long term evolution of the QPO frequencies in 
the atoll source 4U 1 636-53 over an eighteen month period and on the results of their previous research. The 
authors now conclude that in fact there are no peaks in the frequency distribution of the lower kHz QPO in 
this source. In keeping with their previous argument, they conclude that there are no peaks in the frequency 
ratio distribution either. 

While Abramowicz et al. (2003) examine the frequency ratios in pairs of observed frequencies, both 
of the cited papers of Belloni et al. study primarily distributions of the frequencies and focus mainly on the 
lower QPO. Note that in most observations of Belloni et al. (2007) only single QPO frequencies have been 
detected. 

2 Ratio vs. frequency distribution in the atoll source 
4U 1636-53 

Using 4U 1636-53 data from the analysis of Barret et al. (2005b), we prepare a histogram of detections of 
the lower QPO over nine years (from 1996 till 2005) of monitoring by RXTE (Figure 1, left), as well as 
of the upper QPO (Figure 1, right). We restrict our study to frequencies > 400 Hz for the lower QPO, and 
> 800 Hz for the upper QPO. The inaccuracies caused by the long-term frequency drift inside of continuous 
data segments are not important for the purposes of our paper. 

The data of Barret et al. (2005b) have been obtained through a shift-add procedure carried out on indi- 
vidual continuous segments of observation. In this approach each continuous data segment (corresponding 
to few tens minutes of an effective subset from 1 .5 hour RXTE orbital period) is divided into N intervals, and 
searched for a QPO. The shortest usable integration time is estimated such that the QPO is detected above 
a certain significance in at least 80% of the N intervals; a linear interpolation is used to estimate the QPO 
frequency in the remaining intervals. The N PDS are frequency- shifted to the mean QPO frequency over the 
segment and averaged (Mendez et al. 1998). The resulting averaged PDS representing the complete continu- 
ous segment is then fitted with one or two Lorentzians plus a constant corresponding to the counting-statistics 
noise level (Barret et al. 2005b). 

When only one significant peak is detected, the QPO is identified as upper or lower from the parameters 
of the Lorentzian and we refer to such peaks as single QPOs. We stress that the value of the QPO frequency 
itself is not used to distinguish between the upper and lower QPOs — in principle, a QPO of a given frequency 
could be either the upper or the lower QPO. For instance, the quality factor for the lower kHz QPO is a well 
determined function of the frequency, and a different function of frequency for the upper QPO, and we can 
use this and other relationships to identify a QPO of given frequency (see Barret et al. 2005a,b,c, 2006, for 
details). 1 Of course, when two significant kHz QPOs are detected, the upper QPO is the one with the larger 



This way of QPO identification differs from the method based on hardness diagram applied in Belloni 
et al. (2005, 2007). The two methods have a different range of applicability but give comparable results (see 
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frequency, by definition. 

Consequently, the frequency values are averaged through intervals of predetermined length ~ 2000 s. 
Belloni et al. (2005, 2007) analyzed segments of different lengths, resulting in a larger number of detections. 
Hence the histograms we use here are not comparable in details to those of Belloni et al., even for the same 
RXTE data. 

We take into account only the detections of oscillations with quality factor (defined as the QPO centroid 
frequency over the full-width of the peak at its half-maximum) Q > 2 and significance (defined as the 
integral of the Lorentzian fitting the peak in PDS divided by its error) S > 3. 

2.1 Different distributions 

Each of the histograms in the Figure 1 clearly reveals an accumulation of frequencies in the v ~ 900 Hz 
vicinity of the power spectrum. However, for the lower kHz QPO this range (v L ~ 900Hz) is in the high- 
frequency part of the frequency distribution of this QPO, while most detections of Vu are accumulated in the 
low-frequency part of its own range of variation. 

If one is interested (for whatever reason) in the distribution of the frequency ratio v v /v L , then those 
observations in which both QPO peaks are simultaneously detected should be considered. Accordingly, we 
apply our selection criteria to simultaneous significant detections of both QPO frequencies as well. The 
histograms of the upper QPO frequency in this sample (darker bars in Figure 1, right panel) are strikingly 
different from the previous histogram of significant detections of the upper QPO (bars of lighter shade in the 
same figure). A new cluster of frequencies appears, in the range ~ 1 100 Hz to ~ 1200 Hz, at the expense of 
frequencies below 900 Hz, whose occurrence is greatly diminished. 

While there is a positive correlation between the QPO frequencies (e.g., Abramowicz et al. 2005; see 
also Belloni et al. 2005, Zhang et al. 2006), 

vu «0.7v L + 520Hz, (1) 

very clearly the examined data do not support the assumption of Belloni et al. (2005, 2007) that the dis- 
tribution of the ratio of two linearly correlated frequencies is determined by the distribution of one of the 
frequencies even when the second frequency remains undetected — there is apparently no direct link between 
the histogram of all the lower QPO detections (Figure 1, Jeff; lighter) and the histogram of the same QPO 
taken from the subset of twin peak QPO detections (the same figure; darker). This result should have an 
impact on the theory of QPOs. Although a full discussion is beyond the scope of this paper, we note that 
the change in the frequency distribution when a second QPO is detected may be suggestive of a physical 
mechanism, such as mode-coupling. 

To quantify this effect we plot the cumulative distributions of the lower and of the upper QPO, which are 
shown in the Figure 2. Using the Kolmogorov-Smirnov (K-S) test we compare the frequency distributions 
of each (the upper and the lower) QPO measured in all detections, with those measured for the same QPO 
when both the upper and the lower QPO are detected. We obtained the K-S probabilities Pl.ks = 2.35 X 10~ 5 
and pu,KS = 2.24 X 10~ 3 in the case of the lower and upper QPO respectively. Indeed, the two distributions 
are different in both cases. We directly conclude that detections of individual QPOs alone cannot be used for 
calculation of the distribution of the frequency ratios. 

It is interesting to note that the single upper QPOs are mostly detected at relatively low frequencies, 
while the single lower QPOs are detected at relatively high frequencies. Taking into account the linear 
correlation among QPO frequencies, the distributions of single lower and upper QPOs appear to be comple- 
mentary, in the following sense. The lowest-frequency detection of the single lower QPO is at v L = 651 Hz, 
which in the linear correlation corresponds to v D = 976 Hz while the highest-frequency detection of the 
single upper QPO is at v v = 961 Hz, which corresponds to y L = 628 Hz. In other words, if one assumed that 
each of the single upper QPOs is accompanied by a lower QPO of frequency determined from eq. (1), the 
resulting points would all fall to the left of the 3:2 line in Fig. 4 (left panel), and if the same procedure were 
applied to the single lower QPOs, the resulting points would fall to the right of the 3:2 line in the same figure. 
This is illustrated in Fig. 3. Given this fact, it is not surprising that the distribution of actually detected upper 
(or lower) kHz QPOs is completely different from the distribution that would be predicted on eq. (1) from 
the distribution of the other kHz QPO, when detections of single QPOs dominate the data set (Fig. 2). 



e.g., Barret et al. 2005b,c). 
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Figure 2: The cumulative distributions for the kHz QPOs corresponding to Fig. 1, Left: Solid curves 
correspond to the detected lower QPOs (see left panel of Fig. 1), the dotted line labeled "inferred lower" 
indicates the lower QPO frequency calculated from eq. (1) using all detections of the upper QPO. The dashed 
vertical line shows the greatest difference Z) max = 0.515 between the distributions of "all lower" and "lower 
in twin". Right: Analogous lines for the upper QPO. The dashed vertical line on right panel corresponds to 
the maximal difference D mdx = 0.4364 between the "all upper" and "upper in twin" distributions. 
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Figure 3: Distribution of single QPOs. Frequency axes are aligned according to the correlation of eq. (1). 
The shadow denotes a 50 Hz scatter about the lower QPO frequency of 650 Hz, corresponding to a 3:2 ratio. 



2.2 Possible peaks in the ratio distribution 

The left panel of Figure 4 depicts the mutual dependence of frequencies of the lower and upper QPO when 
both were significantly detected. It also display a corresponding histogram of the frequency ratio. As for 
Sco X-l (Abramowicz et al. 2003), this histogram is peaked close to the 3/2 value, and is suggestive of the 
existence of a second peak. 

We have fitted the distribution of the frequency ratios by the sum of two suitably normalized Lorentzians, 

P2(r) = f- o +(!-/) 7 T5 7J' ^ 

where r = v u /v L is the frequency ratio and r\, f"2, A\, A 2 and / are free parameters. Their values obtained 
by the maximum likelihood method are r\ = 1.52, r 2 = 1.28, A { = 0.0327, A 2 = 0.0913 and /' = 0.722, 
reaching K-S probability P2,KS = 0.918. The best fit by a single Lorentzian, with Tq = 1.50 and Aq = 0.0597 
(dotted line in Figure 4) gives the K-S probability Pi.ks = 0.340. Both fits are acceptable. In the right panel 
of Figure 4 we compare cumulative distributions of the observed frequencies with both double and single 
Lorentzians. 
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Figure 4: Left: The frequencies of detected twin QPOs. The inset shows a corresponding histogram of the 
frequency ratio. Right: Cumulative distribution of the frequency ratios, the thick solid line denotes the best 
fit by a sum of two Lorentzians. Fit by a single Lorenzian is marked by dotted line. 



3 Conclusions 

We have demonstrated for a set of uniform data (Barret et al. 2005b) that the frequency distribution of a 
single kHz QPO is not equivalent to the distribution of the corresponding frequency when a pair of kHz 
QPOs have been detected. 

We stress that if there is a one-to-one correspondence between the frequencies and their ratio, as is the 
case for linear functions with a non- vanishing intercept, the question whether to consider the QPO frequency 
distribution or the ratio distribution as fundamental is one of theoretical assumptions, as the two distributions 
are mathematically equivalent. However, the distribution of a single kHz QPO frequency is not predictive 
of the distribution of two frequencies detected simultaneously, nor of the distribution of their ratio, even if 
these frequencies are correlated when both are actually detected. Thus, the study of Belloni et al. (2007), 
who conclude that "there is no preferred frequency or frequency ratio in 4U 1636—53" is based on an invalid 
assumption, and cannot be accepted as applying to the distribution of ratios, as long as it is based on the 
detection of a single frequency. 

The finding that the frequency distribution of a QPO depends on whether or not a second QPO can 
be detected as well should restrict models of the physical origin of the QPO and of X-ray flux modulation, 
regardless of whether or not the value of the frequency ratio is clustered about the specific value of 3/2. 
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